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Summary 


^"■'An  inverse  method  was  developed  for  treating  gas-particle  supersonic  flow 
past  axisymmetric  blunt  bodies.  This  method  is  based  on  two  transformations 
(von  Mises  and  an  additional  one),  which  are  convenient  for  determining  the 
shock-layer  flow  fields  and  the  body  shapes. 

In  using  the  present  method,  the  pure  gas  flow  fields  around  spheres  were 
first  solved  numerically  for  the  freestream  Mach  numbers  h  =  10,  6,  4,  3,  2 
and  1.5.  These  were  found  to  be  in  very  good  agreement  witn^the  available 
results  of  Van  Dyke  and  Gordon.  Then  the  gas-solid-particle  flow  in  the  shock 
layer  around  blunt  bodies  (nearly  spheres)  were  solved  for  the  freestream  Mach 
numbers  10  and  1.5,  with  freestream  loading  ratios  0,  0.2,  0.5  and 

1.0  and  particle  diameters  dpQ=|  1,  2,  5  and  10  uny,  respect ively . .Hie  effects 
of  M^,  dp  and  a  on  the  shock-layer  thickness  ana  the  body  surfacK~pressures  are 
discussed.  The  variations  of  the  flow  properties  along  the  stapiiltion  and 
adjacent  streamlines  are  also  shown  in  some  detail.  \ 
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Notation 


1/2 

freest  ream  gas  sound  speed  =  (iRT  )  , m/sec 

bluntness  of  body  ..hape 
bluntness  of  shock-ware  shape 
drag  coefficient  of  a  particle 

gas  specific  heat  at  constant  pressure,  C  /C  =1 

P  P 

gas  specific  heat  at  constant  volume,  C^/C  =  1/y 

specific  heat  of  solid-particle  material,  ^/Cp  =  fi/'i 

dimensional  particle  diameter,  urn 

particle  diameter,  d  /R 
p  s 

dimensional  gas  thermal  conductivity,  J/(m-sec-K) 

gas  thermal  conductivity,  k/(p  0  C  R  ) 

-  o>  oo  p  s 

-  j 

mass  of  a  particle  =  ^--/h-dp  ,  g 

equilibrium  freestream  Mach  number  [see  Eq.  (5. Ml] 

freest  ream  Mach  number.  0  /a 

number  of  particles  in  a  unit  volume.  1/cm3 


Nusselt  number  of  a  particle 
gas  pressure,  p/c  Uj" 

-  -  2 

stagnation  pressure,  p  /n^U^, 

freestream  gas  pressure,  kPa 

gas  Prandtl  number 

2  2 

gas  constant,  m  /sec  K 

nose  radius  of  body  shape,  R. /R 

dimensional  nose  radius  of  shock  wave  shape,  cm 

nose  radius  of  shock  wave  shape,  R^/R^  =  1 

particle  Reynolds  number,  pwd  /u 

gas  temperature,  T/ (U^ /C  ) 

particle  temperature,  T  / (0m  /Cp) 

freestream  gas  temperature,  K 

gas  velocity  component  in  x-direction,  u/U^ 

particle  velocity  component  in  x-direction,  Up/U^. 

freestream  gas  velocity  =  M  -ao,  m/sec 


gas  velocity  component  in  r-direction,  v/U 


V 

p 

V 

V 

P 

w 


x,  r 


a 


e 

Y 


,5 


in 


r,  '1 


U 


“P 
1  s 


I K  Y 


*P 


V 


P 


particle  velocity  component  in  r-direction  v  /U 

’  p  oo 

transformed  gas  velocity,  V  =  v/Y 

transformed  particle  velocity,  V  =  v  /Y 

P  P 

relative  velocity  between  gas  and  particle 
[(u-up)2  *  (v-vp) 2] 1/2 

cylindrical  coordinates  (see  Fig.  1)  x  =  x/R 
r  =  r/Rs  ’  s’ 

loading  ratio  in  freestream  =  ratio  of  particle 
phase  density  to  gas  density  in  freestream,  p  /p 

p  00 

specific  heat  ratio  of  two  phases,  r  /C 

s  v 

gas  specific  heat  ratio,  C  /C 

P  v 


specific  heat  ratio  of  equilibrium  gas-particle 
mixture  (see  Eq.  (5.13)] 

impact  angle  (see  Fig.  10(a)] 

density  ratio  across  shock  wave  (see  Eq.  (5.11)] 
shock  stand-off  distance,  i/R 

s 

integration  step  size  in  n-direction 
step  size  in  t’ -direction 
coordinates  defined  by  Eqs.  (4.3)  and  (4.4) 
gas  viscosity,  poise  =  g/(cm-sec) 
gas  density,  p/p 

oo 

dimensional  density  of  particle  phase  =  n-m 
density  of  solid  particle  material,  g/cm3 


freestream  gas  density  =  p^/RT^,  g/cm3 
coordinates  defined  by  Eqs.  (3.5)  and  (3.6) 
gas  phase  stream  function 
particle  phase  stream  function 


Superscript 

(  )  dimensional  quantity 


Subscripts 
b  body 

P  particle 

s  shock  or  solid  particle 

st  stagnation 

freestream 


vi 


/ 


1 .  INTRODUCTION 

High-speed  gas  flows  including  small  solid  or 
liquid  particles  have  been  studied  by  a  number  of 
authors  [Refs.  1-b]  in  the  last  few  years  owing  to 
important  engineering  applications  relating  to 
rocket-nozzle  flow,  shock  waves  in  a  dusty  gas 
[Refs.  7,  8]  and  in  supersonic  flight  [Refs.  9-13). 
In  such  two-phase  flows  the  gas-phase  flow  field 
will  differ  from  its  corresponding  field  in  a  pure 
gas  owing  to  the  gas -particle  interaction  through 
viscous  drag  and  heat  transfer. 

The  present  report  deals  with  supersonic  gas 
particle  flow  past  blunt  bodies.  Assuming  that 
the  mass  ratio  of  the  particle  phase  to  the  gas 
phase  in  the  freestream,  or  the  loading  ratio  in 
the  freestream  a,  is  small  enough  and  that  the 
particle  phase  does  not  affect  the  gas  phase, 
I’robstein  and  Fassio  [Ref.  9| ,  Waldman  and  Reinecke 
[Ref.  10),  and  Spurk  and  Gerber  [Ref.  11)  calculated 
the  particle  trajectories  in  the  hypersonic  shock 
layer.  Their  main  interest  was  to  find  the  collec- 
■  ■/’:  ■  •  \  v,  defined  as  the  ratio  of  the  number 
of  particles  which  strike  the  body  to  the  number 
of  particles  which  could  strike  it  if  their  trajec¬ 
tories  were  straight  and  unaltered  by  the  vehicle 
flow  field. 

Taking  into  account  the  interactions,  i.c.,  the 
drag  force  and  the  heat  exchange  between  the  gas 
and  the  particles,  Chang  [Ref.  13)  investigated 
dusty  supersonic  flow  past  blunt  bodies  using  a 
direct  method.  He  showed  that  the  shock-layer  flow 
is  in  nonequilibrium  and  investigated  the  effects 
of  the  loading  ratio  a  on  this  flow.  His  main 
interest  was  to  clarify  the  nonequilibrium  pheno¬ 
mena  in  the  two-phase  flow  within  the  shock  layer, 
which  is  also  the  purpose  of  the  present  report. 

In  this  report  a  numerical  method  for  treating 
dusty  supersonic  flow  past  blunt  bodies  will  be 
developed.  This  method  is  the  so-called  inverse 
method  in  which  the  shock  shape  is  specified  and 
the  flow  field  and  body  shape  is  to  be  determined. 
First,  the  governing  equations  for  dusty  flow  past 
axisymmetric  bodies  are  derived  in  a  cylindrical 
coordinate  system  (x,  r) .  It  is  assumed  that  there 
are  no  collisions  between  the  particles.  The  inter¬ 
actions  between  the  gas  and  the  particles  considered 
are  only  the  drag  force  and  the  heat  exchange. 

Next,  the  governing  equations  are  transformed  into 
a  (*,  Y)  coordinate  system  by  using  a  von  Mises 
transformation,  where  ti  is  the  stream  function 
for  the  gas  phase,  and  Y  =  r.  To  facilitate  the 
integration  of  the  governing  equations  in  the  shock 
layer,  they  are  further  transformed  into  a  (£,,  n) 
coordinate  system,  which  will  be  discussed  subse¬ 
quently.  The  advantages  of  the  present  method  are: 

(1)  body  shapes  are  easily  determined  without  any 
approximate  calculations;  (2)  relatively  low  com¬ 
puter  storage  is  required,  since  it  is  unnecessary 
to  store  all  the  flow  field  data  within  the  shock 
layer;  and  (3)  solutions  are  obtained  quickly. 

Although  the  advantages  and  disadvantages  of 
the  inverse  method  and  the  direct  method  (including 
the  time-dependent  method)  for  the  blunt-body 
problem  have  been  fully  discussed  in  Refs.  14-19, 
it  is  worthwhile  to  summarize.  The  direct  method 
is  more  important  for  practical-type  problems  than 
the  inverse  method.  However,  the  inverse  method 
is  significant  and  useful  for  complex  flow  problems 


which  have  not  been  fully  clarified,  such  as  three- 
dimensional  flows  of  a  pure  gas,  and  two  or  three- 
dimensional  flows  with  relaxation  processes.  The 
inverse  method  is  much  easier  for  numerical  calcula¬ 
tions,  gives  a  more  detailed  flow  field,  and  provides 
better  insight  into  the  flow  fields. 

Using  the  present  method,  the  dusty-gas  shock 
layer  flows  around  hHmt  bodies  (nearly  spneics) 
were  solved  numerically.  For  experimental  purposes, 
glass  microspheres  in  air  were  chosen  as  the  dusty- 
gas.  To  investigate  the  effects  of  the  freestream 
Mach  number  M„,  particle  diameter  dp,  and  freestream 
loading-ratio  a  on  the  shock  layer  of  the  dusty-gas 
flow,  the  following  parameters  were  chosen: 

M^  =  10,  1.5 

d  =  1 ,  2 ,  5 ,  1 0  urn 
P 

a  =  0,  0.2,  0.5,  1.0 

Numerical  results  for  a  pure  gas  (a  =  0)  were  com¬ 
pared  with  those  of  Van  Dyke  and  Gordon  [Ref.  lb] 
and  very  good  agreement  was  obtained,  thus  lending 
confidence  to  the  dusty-gas  calculations. 

It  is  worth  noting  that  spalling,  ablation  and 
erosion  of  the  glass  particles  were  not  considered 
in  the  present  study. 

2.  GOVERNING  HQUAT10NS 

To  formulate  the  motion  of  the  gas-particle  flow 
in  the  shock  layer,  the  following  assumptions  are 
made  [Ref.  4): 

(1)  The  volume  fraction  of  the  particles  is  small  so 
that  the  collisions  between  individual  particles 
are  neglected. 

(2)  The  interactions  between  the  gas  and  particles 
consist  only  of  the  drag  force  between  the  gas 
and  particles  and  the  heat  exchange  between 
them. 

(3)  The  particles  are  assumed  to  be  spheres  of  uni¬ 
form  size. 

(4)  The  number  of  particles  is  large  and  the  flow 
can  be  treated  as  a  continuum. 

In  this  report  a  superscript  (  )  over  a  quantity 
means  a  dimensional  quantity,  otherwise  it  is  non- 
dimensional.  Let  u  and  v  be  the  x  and  r  components 
of  velocity  (Fig.  1),  o  the  density,  p  the  pressure, 
and  T  the  temperature  of  the  gas,  and  up,  v_,  Pp 
and  Tp  be  the  corresponding  values  of  theparticle 
phase.  Here,  pp  =  Pp/p^  =  n-ra/P^,  where  n  is  the 
number  of  particles  per  unit  volume  and  m  is  the 
mass  of  a  particle.  All  lengths  are  nondimensional - 
ized  by  the  shock  nose  radius  Rs,  velocities  by  the 
freestream  velocity  U^,  densities  by  the  freestream 
gas  density  o  ,  pressure  by  p  UB‘,  and  temperatures 
by  U„2/c  where  Cp  is  the  gas  specific  heat  at 
constant  pressure .  The  other  nondimensional  quan¬ 
tities  are  listed  in  the  Notation. 


Under  the  above  assumptions,  the  governing 
equations  of  the  steady,  axisymmetric,  dusty-gas 
flow  become  (Refs.  4,  13): 
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For  the  gas  phase. 
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Momentum : 
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(2.1) 
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'  f  u 
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J  P  px  P  pr 

(2.11) 

that  is,  the  kinetic  energy  of  the  particles  is  only 
affected  by  the  work  done  by  the  drag  on  the  parti¬ 
cles. 


The  drag  force  exerted  by  n-part  icles  of  dia¬ 
meter  dp  in  a  unit  volume  (=RSJ  )  of  space  moving 
through  a  perfect  gas  is  given  by 
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Fnergv : 
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State: 
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and  for  the  particle  phase, 


Font inuitv : 


is  the  relative  velocity  between  the  gas  and  part 
cles,  and  Cp  is  the  drag  coefficient.  At  present 
(2.5)  the  nonstationary  drag  coefficient  for  a  gas-part] 
cle  flow  is  uncertain  (Ref.  6],  and  is  assumed  as 
(Ref.  20] 


CD  =  0.48  ♦  28/Re 


0.85 


(2 


where 


—  (r.  u  )  ♦  —  (r.  v  )  =  0  (2.6) 
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(2.16) 


Momentum: 
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Energy : 
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where  cs  is  the  specific  heat  of  the  particles. 

FpX  and  Fpr  are  the  components  of  the  force  exerted 
upon  a  unit  volume  of  gas  by  the  particles,  and  Qp 
is  the  heat  transferred  from  the  particles.  The 
terms  UpFpX  +  vpFpr  represent  the  work  done  by  the 
particles  passing  through  the  gas.  Equation  (2.4) 
shows  that  the  thermal  energy  of  the  gas  is  changed 
by  both  the  heat  exchange  and  the  dissipation  due 
to  the  particles  moving  relative  to  the  gas.  Equa¬ 
tion  (2.9)  states  that  the  thermal  energy  of  the 
particles  can  be  changed  only  by  heat  transfer. 

From  Eqs.  (2.2)  and  (2.3), 


2  2 
u  ♦  v 


-u  -  v  ^  *  uF  ♦  vF 
7x  Tr  px  pr 


is  the  Reynolds  number  based  on  the  diameter  of 
the  particle  dp  and  the  relative  velocity  w. 

The  heat  transfer  between  the  gas  and  particles 
is  given  by 


Q  =  n  Nu  ■  d  k(T  -T) 
p  p  p 


(2.17) 

where  the  Nusselt  number  is  expressed  as  [Ref.  21] 

(2.18) 


Nu  =  2.0  ♦  0.6  Pr1/3  Re1/2 


(2.19) 


where  Pr  is  the  Prandtl  number, 
uc 

Pr  = 

k 

and  k  is  the  gas  thermal  conductivity. 


In  the  present  analysis,  the  Prandtl  number  is 
taken  at  0.75,  the  specific  heat  of  the  gas  is  held 
constant,  the  specific  heat  ratio  of  the  two  phases 
is  taken  as  8  =  cs/cv  -  1.0,  and  the  power  law  for 
the  viscositv-temperature  variation  is  assumed  as 
[Ref.  22] 


(2.10) 


U  =  1.71  x  10 


-4 


T 

273 


0.77 


poise 


(2.20) 


Similarly,  from  Eqs.  (2.7)  and  (2.8), 


It  is  worth  noting  that  any  other  relation  for 


2 


/ 


viscosity  can  be  easily  inserted  without  affecting 
other  parts  of  the  analysis. 


p  =  —•  PT 


3 .  EQUATION  OF  MOTION  IN  (y,  Y)  COORDINATE  SYSTEM 


Similarly,  the  governing  equations  (2.b)-(2.9)  for 
the  particle  phase  arc  transformed  to 


Introduce  the  stream  functions  ;>(x,  r)  for  the 
gas  phase  and  ,p(x,  r)  for  the  particle  phase  as 


'll  -V  'V  , 

p  P  pi  p 

V  — A-  ♦  -Tif-  +  .  — L  *  —  lu  ,  *  v  .  : 

x  ",  A  r  p  x  p  r 

P  1 

v  •.  v 

♦  ^  -  o 


11  'll 

.  ( u  ■_  +  v  .  )  — 4  +  ■  v  — =  -  T  I  3  .  1 8  j 

P  p  x  p  r  P  P  A  px 


•V  'V 

.  (u  .  ♦  v  .  )  —1 2  ♦  .  V  -rrP  =  -1;  (3.191 

p  p  x  p  r  •.  p  p  A  pr 


which  satisfy  the  continuity  equations  (2.1!  and 
(2.0).  Using  the  von  Mises  transformation  (Ref. 

23]  to  transform  the  governing  equations  in  (x,  r) 
coordinate  system  into  those  in  a  (,,  Y)  coordinate 
system,  where  Y  =  r,  a  set  of  transformation  oper¬ 
ators  can  he  obtained  as  [Refs.  24,  25 ] 


•T  -T 

.  (u  .  ♦  v  .  1C.  — 4  *  .  v  C  —tP  =  -Q  (3.20) 

p  p  x  p  r  s  ■.  p  p  s  A  p 


4  .  EQUATIONS  OF  MOTION  IN  ( ■'  ,  •)  COORD  1  NATL  SYSTEM 

The  shock-layer  gas-particle  flow  is  now  analysed 
by  the  inverse  method,  that  is,  the  shock  wave  is 
prescribed  and  the  body  surface  is  then  determined. 


•r  -Y  'r  <•;. 


Tor  the  stream  functions  ahead  of  the  shock  wave, 
choose 

..4 


From  the  inverse  transformation, 


•Y  r  A'  ex 


To  facilitate  the  integration  of  the  governing  equa¬ 
tions  in  the  shock  layer,  introduce  further  the 
following  independen;  variables  [Ref.  24]: 


Comparing  Eqs.  (3.")  and  (3.81  and  Eqs.  (3.5)  and 
(3.9)  , 


Using  the  above  transformation  operators,  the 
governing  equations  (2.1)-(2.5)  for  the  gas  phase 
are  transformed  as  follows: 


In  this  ( : ,  *,)  coordinate  system,  the  shock  wave 
surface  and  the  body  surface  correspond  to  h  =  1  and 
0,  respectively.  From  the  transformation  equations 
(4.3)  and  (4.4),  a  set  of  transformation  operators 
are  obtained: 


Y  AY  A£  ?,  3r 


VC  -TT  -  V 


'V  Tv  V  Ap  V  „ 

+  Vrly*vT*p^Y*Y  =  0 

(3.12) 
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1  A 

t  An 

(4.b) 

'  v  44  =  -v  .  F 

•Y  x  'll  px 

(3.13) 

Considering  v  «■  Y,  Vp 
dependent  variables  v 

"  Y  near  Y  =  0, 
and  Vp  to  V  and 

transform  the 
Vp  as 

■  44  =  -  -  y  4£  .  f 

<Y  Y  r  ?’i‘  pr 

(3.14) 

V 

v 

Y 

(4.7) 

■  44  =  Q  *  (u  -u)F  ♦  (v  -v)F 
>1  Mp  p  px  P  Pr 

(3.15) 
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Y 

(4.8) 
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lhon  the  governing  equations  (3. 1.2) -(3. 20)  become, 
for  the  gas  phase. 


•u  •*  ;V 


■'  !  1  -v  1  ...  \  1 


=  Jii  _  ii 
,v  >oV 


(■).;») 

(4.10) 


5 .  IN1T1AI.  CONDITIONS 

The  initial  conditions  immediately  behind  the 
shock  wave  for  the  gas  are  the  usual  jump  conditions 
( Rank  ine-llugon  iot  relations),  and  for  the  dust 
particles  use  the  frozen  eonditions,  that  is,  the 
particle  flow  quantities  immediately  behind  the 
shock  wave  are  the  same  as  those  immediately  ahead 
of  the  shock  wave. 
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he  the  equation  of  the  shock -wave  shape.  Then  the 
flow  quantities  behind  the  shock  wave  are  given  by 
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and  for  the  particle  phase, 
U' 
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for  the  gas,  and 
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liquation  (3.11)  becomes 
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doV 


(4.18) 


u  =  1 
P 


v  =  n 

p 


'p 


for  the  particles,  where 


i-i  .  s" 
1  "  i  +  l  w  „d 


In  deriving  the  above  equations,  the  following  rela¬ 
tions  were  used: 


1  *  s 


(5.') 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.1.’) 


-T-  =  -2,'V 


(4.19) 


T  V-Y‘n 


=  d(ou-n)  =  4,V\ 


■>x 


(4.20) 


For  the  ten  unknown  variables,  u,  V,  .  ,  T, 
p,  Up,  Vp,  0p,  Tp  and  x,  the  ten  equations  (4.9)- 
(4.18)  provide  a  solution  through  numerical 
integrat ion . 


and  y  Is  the  specific  heat  ratio  of  the  gas,  and 
M,  »  U,„/a,„  is  the  freestream  Mach  number. 

It  is  worth  noting  that  the  specific  heat  ratio 
of  the  equilibrium  gas-particle  mixture  >c  and  the 
equilibrium  Mach  number  ahead  of  the  shock  wave  Mc 
are  given  by  (Ref.  b] 


‘■(+0*1. 

1  ♦.»  •  h 


(5.13) 
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where  i  is  the  loading  ratio,  and  B  =  cs/cv  is  the 
ratio  of  the  specific  heats  of  the  two  phases. 


b.  NUMERICAL  PROCEDURE 


The  governing  equations  (4 .9) -(4 . 18)  are  numer¬ 
ically  integrated  from  a  point  immediately  behind 
the  shock  wave  (n  =  1)  to  the  body  surface  (n  -  0) 
by  a  forward  (i.e.,  downstream)  integration  method 
[Ref.  16]. 

The  numerical  procedure  is  summarized  as  follows: 

1.  Calculate  the  initial  values  f ' m '  at  '  =  1  from 
F.qs .  IS  .2)  —  (5.10),  where  f(m)  means  any  flow 
quantity . 

2.  Calculate  •'  derivatives  f.-m^  by  numerical 
differentiation  (here  a  S-point  Lagrange 
relation  was  used) . 

3.  Calculate  f^  from  the  governing  equations 
(4.9)  —  (4. 18). 

4.  Extrapolate  linearly  to  the  next  smaller 

value  of  ri : 

f(m*D  =  f(m)  .  in.f  (m)  (f> . , , 


5.  Calculate  f,  derivatives  f,-m+^  as  jn  step  2. 
b.  Calculate  f^1"*1-*  as  in  step  3. 

7.  Extrapolate  again  f(m+1)  using  average  value  of 
f  n  * 

f(m+l)  _  (m)  _  An  ff(m)  +  f(m*l) 

2  In  n 

8.  Repeat  steps  2  to  7  for  this  new  value  of  r, 
and  continue  until  n  =  0. 

As  the  freestream  Mach  number  M  reaches  a  low 
value,  the  numerical  calculations  become  difficult 
owing  to  numerical  instability  at  the  outer  calcul¬ 
ation  points  (end  instability).  In  this  case,  the 
end  stability  is  eliminated  by  dropping  the  outer¬ 
most  2  or  3  points  at  each  step. 

The  mesh  sizes  AF),  An  used  were  AF,  =  0.005,  0.01 
and  An  =  0.1,  0.2.  The  numerical  calculations  were 
done  on  an  IBM  1130  computer  at  UTIAS,  and  the 
calculation  time  was  about  1  to  3  seconds  for  each 
case . 


(6.2) 


7.  NUMERICAL  RESULTS  AND  DISCUSSIONS 
7.1  Pure-Gas  Case 


To  check  the  adequacy  of  the  present  numerical 
results,  the  flow  field  around  a  sphere  in  a  pure 
gas  (i.e.,  no  solid  particles)  was  solved  and  com¬ 
pared  with  those  of  Van  Dyke  and  Gordon  [Ref.  16], 

Any  shock  wave  ahead  of  a  smooth,  blunt  body 
may  be  described  by  [Refs.  15,  16] 


l 

j 

(7.1a)  j 
or 

r“  =  2x  -  B  x“  (7.1b) 

s 

here  Rs  is  the  nose  radius,  Bs  is  the  bluntness  of 
the  shock  shape,  and  r  =  r/Rs,  x  =  x/Rs.  The  blunt  - 
ness  Bs  is  a  convenient  parameter  that  characterizes 
the  eccentricity  of  the  conic  section.  It  is  worth 
noting  that  Bs  •  0  generates  a  hyperbola,  Bs  =  0 
generates  a  parabola,  and  Bs  0  generates  an  ellipse, 
and  Bs  =  1  results  in  the  special  case  of  a  sphere. 
Table  1  shows  the  values  of  Bs  taken  from  Ref.  16. 

In  addition,  other  quantities  such  as  shock  standoff 
distance  A,  nose  radius  Rb  and  bluntness  Bb  are  also 
tabulated  and  will  be  discussed  subsequently. 


=  2R  x 


B  x 
s 


The  body  shape  corresponding  to  a  given  shock 
shape  may  be  approximated  by  [Refs,  lb-19] 

r~  =  2Rj;(x-A)  -  Bfa(x-A)‘  (7.2a) 

or  ,  , 

r"  =  2Rb(x-.’.)  -  Bb(x-.‘,r  (7.2b)  i 

where  Rb  =  Rh/Rs  is  the  nondimensional  nose  radius, 

Bb  the  bluntness  of  the  body,  and  A  -  A/Rs  the  non- 
dimensional  standoff  distance.  The  process  used  to 
determine  Rb  and  Bb  is  the  same  as  that  of  Refs,  lb 
and  17,  that  is,  the  quantity 


2  (x -A) 


is  plotted  versus  <x-A),  and  the  points  are  approx¬ 
imated  by  a  straight  line.  Then,  the  vertical 
intercept  of  the  line  is  Rb,  and  its  slope  is  -[Bb- 

Figures  2(a) -(d)  show  the  shock  waves,  sonic 
lines  and  streamlines  for  a  sphere  at  freestream 
Mach  numbers  Mr  =  10,  3,  2  and  1.5,  respectively. 

The  abscissa  x  =  x/Rs  and  the  ordinate  r  =  r/Rs  are 
nondimensional  distances,  as  noted  previously.  For 
comparison,  Van  Dyke  and  Gordon's  numerical  results 
[Ref.  lb)  are  shown.  The  present  calculated  body 
shapes  agree  very  well  with  those  of  Ref.  lb  for  the 
high  freestream  Mach  number  M^  =  10,  and  the  agree¬ 
ment  is  also  good  for  the  low  freestream  Mach 
numbers  M  =  3,  2  and  1.5.  The  present  calculated 
sonic  line  shapes  almost  coincide  with  those  of 
Ref.  lb  except  for  the  case  of  Mai  «  1.5.  It 
may  be  noted  that  the  present  sonic-line  shape  for 
M^  =  1.5  resembles  that  of  Ref.  27.  The  flow 
patterns  of  the  shock  layer  around  a  sphere  are 
well-defined  from  the  present  calculated  stream- 
1 ines . 

In  Figs.  3(a)  and  (b)  appear  the  relations 
between  the  quantities  r‘/2(x-A)  and  (x-A)  for  the 
calculated  bodies.  From  these  figures  Rb  and  Bb 
are  determined  as  shown  in  Table  1.  It  is  seen  that 
the  present  nose  radii  Rb  coincide  with  those  of 
Van  Dyke  and  Gordon  (Ref.  16]  within  1.5%,  and  the 
present  bluntness  of  the  body  Bb  ~  1,  i.e.,  the 
calculated  body  shapes  are  spheres. 

The  shock  standoff  distances  A  are  also  listed 
in  Table  1.  The  present  results  agree  very  well 
with  those  of  Ref.  16  for  high  freestream  Mach 

numbers  M^  =  10,  6  and  4,  and  also  agree  quite  well 

with  those  of  Ref.  lb  for  low  freestream  Mach 

numbers  M  =  3,  2  and  1.5  within  2.5%. 

oo  ’ 
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Table  2  shows  the  stagnation  pressure  on  the 
bodies.  The  exact  values  shown  in  Table  2  are 
calculated  from  the  following  equation  [Ref.  26]: 


Pst  1  (  l*'1 


'  1 


1  +  1 


2,M“  -  h-1) 


Y 

i 


(7.3) 


The  present  results  coincide  with  the  exact  values 
within  U. 


The  pressure  distributions  on  the  bodies 
[spheres)  are  shown  in  Figs.  4(a)-(d)  for  Mv  =  10, 

3,  2  and  1.5,  respectively.  The  present  results 
agree  very  well  with  those  of  Van  Dyke  and  Gordon 
[Ref.  16]  in  the  stagnation  region  for  all  free- 
stream  Mach  numbers.  But  the  difference  between  the 
present  results  and  those  of  Ref.  lb  appear  near 
sonic  lines  on  the  bodies  at  the  low  freestream 
Mach  numbers  Mr  =  3,  2  and  1.5.  It  may  be  worth 
noting  that  Hamilton  [Ref.  19]  has  pointed  out  that 
the  results  of  Ref.  lb  underpredict  the  pressure 
near  the  sonic  line  on  the  body  at  low  freestream 
Mach  numbers. 


7 . 2  Dusty-Gas  Case 


Eq .  (2.20),  and  the  drag  coefficient,  Eq .  (2.15), 
remained  unchanged. 

Figures  5(a)— (d)  show  the  assumed  shock  waves 
and  calculated  body  shapes  for  Mo  =  10  and  dp  =  1 , 

2,  5  and  10  pm,  respectively.  The  loading  ratios 
i  ahead  of  the  shock  wave  were  »  =  0,  0.2,  0.5  and 
1.0.  The  abscissa  and  ordinate  are  nondimensional - 
iced  by  the  shock  nose  radius  Rs  =  1  cm.  It  is 
seen  that,  when  dp  =  1  am,  the  shock  layer  thickness 
decreases  as  i  increases  from  0  to  0.5  [Fig.  5(a)]. 
(When  a  =  1.0,  a  body  shape  was  not  obtained 
because  of  numerical  instability  near  the  body 
surface.)  This  tendency  is  similar  to  the  results 
of  Ref.  13  for  an  aerosol-gas  mixture.  However, 
when  dp  =  2,  5  and  10  urn,  the  shock  layer  thicknesses 
increase,  that  is,  the  shock  standoff  distance 
increases  with  i.  The  reason  will  be  considered 
later,  since  the  shock  stand-off  distance  is  closely 
associated  with  the  variation  of  the  gas-flow 
quantities  within  the  shock  layer. 

Figures  6(a)-(d)  show  the  relations  between 
the  quantities  r2/2(x-.'.)  and  ( x - )  of  the  calculated 
bodies.  From  these  figures,  the  nose  radius  Rb  and 
the  bluntness  of  the  body  Bfo  were  calculated,  which 
are  shown  in  Tables  4 (a) -(d).  It  is  seen  that,  for 
example,  when  dp  =  1  urn,  Rb  increases  with  increasing 
i,  and  Bb  differs  slightly  from  unity  with  increasing 
i.  It  may  be  noted  that  Bs  =  1.0  means  a  true  sphere, 
and  Bs  =  0.93  means  a  prolate  sphere. 


The  following  numerical  examples  concern  the 
flow  of  gas-solid  particles  past  blunt  bodies 
(near  spheres).  Since  the  inverse  method  is  used, 
the  f reest ream-flow  conditions  and  the  shock  waves 
described  by  Eq.  (7.1)  arc  assumed  first  and  then 
the  numerical  method  is  used  to  find  the  shock-layer 
gas-particle  flows  and  the  body  shapes. 

Air  (i  =  1.4,  gas  constant  R  =  287  m“/s*k)  and 
glass  spheres  (. s  -  2.5  g/cms)  are  used  for  the 
gas-solid  particle  flow.  The  temperature  and  the 
pressure  of  the  air  ahead  of  the  shock  wave  are 
T^  -  300  K,  p^  =  101.3  KPa,  and  the  reference 
length  is  Rs  =  1  cm  for  all  the  calculations. 

To  investigate  the  effects  of  the  freestream 
Mach  number  the  freestream  loading  ratio  a  and 
the  particle  diameter  dp  on  the  shock  layer  gas- 
particle  flow,  the  following  range  of  parameters 
were  used: 


Figures  7 (a) -(d)  show  the  variations  of  the 
velocities  and  the  temperatures  of  the  gas  and 
particles  along  the  stagnation  streamlines  for 
dp  =  1,  2,  5  and  10  urn,  respectively.  Similarly, 
Figs.  8(a)-(d)  show  the  variation  of  the  gas  and 
particle  densities  and  gas  pressures  along  the  stag¬ 
nation  streamlines.  As  noted  previously,  the 
particle  density  ,p  =  Pp/o  *  m-n/f  ,  and  it 
changes  within  the  shock-layer  flow  field. 

The  results  for  dp  =  1  urn  will  now  be  discussed 
in  detail.  Figure  7(a)  shows  that  the  gas  velocity 
u  decreases  gradually  .long  the  stagnation  stream¬ 
line  and  becomes  :ero  at  the  body  surface.  However, 
the  particle  velocity  Up  decreases  much  more  rapidly 
as  a  result  of  the  drag  interaction  of  the  particles 
with  the  gas.  The  particles  reach  almost  the  same 
velocity  as  the  gas  near  the  body  surface.  In  other 
words,  the  particle  and  gas  velocities  are  in 
dynamic  equilibrium  near  the  body  surface. 


M  =  10,  1.5 

l  =  0,  0.2,  0.5,  1.0 

d  =  1,  2,  5,  10  urn 
P 

The  number  of  solid  particles  in  the  freestream 
flow  ahead  of  the  shock  wave  appear  in  Table  3. 


7.2.1  Results  for  M^  =  10 

Although  the  results  for  M  =  10  are  only  of 
academic  interest  because  of  tfie  assumptions  made 
concerning  gas  properties  at  high  temperature  and 
the  resulting  particle  drag  and  heat  transfer 
coefficients,  these  calculations  were  performed  in 
order  to  show  the  trends  of  the  shock-layer  gas- 
particle  flows.  As  noted  previously,  the  Prandtl 
number  was  taken  at  0.75,  the  specific  heat  of  the 
gas  was  held  constant,  and  the  viscosity  power  law. 


Figure  7(a)  also  shows  that  the  gas  temperature 
T  increases  along  a  stagnation  streamline  until  it 
reaches  a  maximum,  and  then  decreases  as  the  loading 
ratio  ahead  of  the  shock  i  increases  from  0.2  to 
1.0.  This  variation  of  the  gas  temperature  is 
explained  as  follows.  As  noted  previously  [Eq.  (2.4)] 
the  gas  temperature  T  is  changed  owing  to  both  the 
work  done  by  the  particles  on  the  gas  (Up-u)FpX  and 
the  heat  transfer  from  the  gas  to  the  particles  Qp. 
Just  behind  the  shock  wave,  the  difference  between 
the  gas  and  particle  velocities  is  so  large  that  the 
heat  generated  by  the  work  done  by  the  particles 
(Up-u)FpX  dominates  the  heat  loss  Qp,  and  as  a  result 
the  gas  temperature  T  increases  behind  the  shock  wave 
at  first.  Near  the  body  surface,  the  velocity  differ¬ 
ence  (Up-u)  becomes  small,  and  (Up-u)FpX  becomes  small 
However,  the  increase  in  Pp  results  in  a  larger  heat 
loss  from  the  gas  to  the  particles.  Consequently, 
the  heat  loss  Qp  dominates  the  work  done  (Up-u)FpX, 
and  the  gas  temperature  T  decreases  near  the  body 
surface.  When  the  heat  generated  by  the  work  done 
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by  the  particle  on  the  pas  is  equal  to  the  heat  loss 
to  the  particle  phase,  the  pas  temperature  reaches 
a  maximum.  This  variation  of  gas  temperature  is 
similar  to  that  obtained  in  Ref.  13. 

It  should  he  noted  that  the  pas  stagnation 
temperature  on  the  body  increases  with  increasing 
i  owing  to  the  decrease  in  kinetic  energy  of  the 
particles.  fhe  particle  temperature  Tp  increases 
monotonical ly  because  of  the  heat  transfer  from 
the  gas  phase,  but  it  does  not  reach  the  gas 
temperature  T.  The  gas  and  particle  temperatures 
are  always  in  nonequilibrium. 

Figure  8(a)  shows  that  the  gas  density  .  for 
i  0  decreases  along  the  stagnation  streamline  at 
first  and  then  it  increases  near  the  body  surface. 
However,  the  density  of  the  particle  phase  .  p 
increases  gradually  at  first  and  then  increases 
more  rapidly  ntar  the  hodv  surface  as  .  increases 
from  0.2  to  1.0.  The  gas  pressure  p  increases 
along  a  stagnation  streamline  owing  to  the  decrease 
of  the  gas  and  particle  velocities. 

The  results  for  dp  =  ..m  [Fig.  8(h)]  are 

significantly  different  than  for  dp  =  1  ,.m.  In 
this  case,  the  particle  velocity  up  does  not  reach 
the  gas  velocity  u  near  the  body  surface  [Fig.  “(h)], 
and  a  substantial  difference  between  the  gas  and 
particle  velocities  Uip-u)  always  exists.  Heat  is 
always  generated  in  the  gas  phase  owing  to  the  work 
done  hv  the  particles.  On  the  other  hand,  the 
particle  temperature  Tp  does  not  increase  very  much. 
That  is,  the  heat  loss  from  the  gas  to  the  parti¬ 
cles  is  slight.  Consequently,  the  gas  temperature 
T  always  increases  along  the  stagnation  streamline. 
When  i  =  0.5  and  1.0,  the  gas  temperature  T  reaches 
more  than  10,000  k  near  the  body  surface.  In  these 
cases,  one  has  to  consider  real-gas  effects,  and 
evaporation  and  melting  of  the  particles.  Conse¬ 
quently,  the  range  of  <  for  which  the  present 
analysis  is  applicable  is  .  0.5. 

Figure  8(h)  shows  that  the  density  of  the  gas 
,  decreases  along  a  stagnation  streamline  owing  to 
the  increase  in  gas  temperature  as  .  increases  from 
0.2  to  1.0.  The  gas  pressure  p  ana  the  density  of 
the  particle  phase  p  increase  along  a  stagnation 
streamline  as  i  increases. 

As  mentioned  previously,  when  dp  -  2,  5  and  10 
cm,  the  shock  standoff  distance  increases  with 
increasing  .  [Figs.  5(b),  (cl,  id)].  This  pheno¬ 
menon  can  now  he  explained  as  follows:  the  gas 
density  .  along  stagnation  and  adjacent  stream¬ 
lines  decreases  owing  to  the  increase  in  gas 
temperature  as  i  increases.  On  the  other  hand, 
the  shock  location  is  determined  by  the  continuity 
condition  of  equal  gas  mass  flow.  Consequently, 
the  shock  standoff  distance  increases  as  i  in¬ 
creases  . 

The  results  for  dp  =  5  and  10  urn  appear  in 
Figs.  7(c|,  (d)  and  8(c),  (d).  As  the  particle 
diameter  dp  increases,  the  rates  of  decrease  of 
Up  and  increase  of  Tp  along  a  stagnation  stream¬ 
line  decrease.  In  other  words,  the  particle  flow 
behaves  almost  like  a  frozen  flow  as  dp  increases. 

In  the  case  of  dp  =  5  and  10  urn,  the  variation  of 
the  other  flow  quantities  u,  T,  p,  ,  and  ,p  along 
a  stagnation  streamline  are  qualitatively  'the 
same  as  for  dp  =  2  qm 


Figures  9(a)-(d)  show  the  pressure  distribu¬ 
tions  on  the  body  surfaces  for  dp  -  1,  2,  5  and 
10  urn,  respectively.  In  each  case,  the  pressure 
increases  with  the  loading  ratio  i,  ahead  of  the 
shock  wave.  However,  for  ■  =  0,  the  agreement 
with  the  pressure  distribution  observed  in  Ref. 
lb  is  very  good. 

Figures  10(a)— (J)  show  the  gas  and  particle 
streamlines  near  the  stagnation  regions  in  the 
shock  layers.  When  dp  =  1  ..m,  the  particle  paths 
are  straight  just  behind  the  shock  wave,  and 
continue  to  deflect  as  they  approach  the  body 
surface  and  collide  at  small  impact  angle  as 
shown  in  Fig.  10(a).  As  dp  increases  from  2  to 
10  ;.m,  the  particle  deflection  decreases  and  they 
tend  to  move  in  a  straight  path  towards  the 
body  surface  where  they  collide  with  large  impact 
angles . 

After  the  particles  collide  with  the  body  sur¬ 
face,  they  may  break  up  or  reflect  or  stick  to  the 
body  surface.  However,  at  present,  it  is  not  known 
how  the  particles  behave  after  colliding.  This 
analysis  may  be  valid  if  the  particles  break  after 
colliding  into  much  smaller  ones  to  form  a  thin 
dusty  layer  or  if  they  stick  to  the  body  surface. 

“.2.2  Results  for  M _  -_  1  . 5 

Figures  Ilia)- Id)  show  the  assumed  shock 
shapes  and  calculated  body  shapes  for  H,  -  1.5 
and  dp  -  1,  2,  5  and  10  ,.m,  respectively.  In  the 
case  of  M.  =  1.5,  the  shock  standoff  distances 
decrease  with  the  loading  ratio  i  ahead  of  the 
shock  wave  for  all  particle  diameters  dp.  The 
greatest  effect  on  the  shock  standoff  distance 
occurs  for  the  smallest  particle  diameter  dp  =  lam. 
lbc  reason  is  that  when  dp  *>  1  i.m,  the  gas  density 
increases  most  significantly  along  the  stagna¬ 
tion  streamline  and  adjacent  streamlines  owing 
to  the  gas  and  particle  interaction  as  i  increases 
| see  Fig.  11(a)).  Consequently,  the  shock  stand¬ 
off  distance  decreases  as  •  increases.  As  the 
particle  diameter  d„  increases  from  2  to  10  cm, 
the  effect  of  .  on  the  shock  standoff  distance 
becomes  weak.  Wien  dp  =  10  urn,  the  effect  of  i  is 
indistinguishable  from  that  of  a  pure -gas  flow. 

•> 

In  Figs.  12(a)-(d),  the  quantities  r“/2(x-.M 
are  plotted  versus  (x-  ) .  In  these  cases,  the 
calculated  points  near  a  sonic  point  deviate 
considerably  from  a  straight  line.  Tables  5(a)-(d) 
show  the  shock  standoff  distance  nose  radius  R(, 
and  hluntness  B|,  of  the  calculated  bodies  for 
dusty-gas  flows.  It  is  seen  that,  for  all  particle 
diamet.  rs  dp,  the  nose  radius  Rp  increases  as  a 
increases  and  the  bluntness  B[,  becomes  smaller 
than  unity,  i.e.,  the  calculated  body  shapes 
become  a  prolate  sphere  as  ■  increases. 

Figures  13(a)-(d)  and  14(a)-(d)  show  the  vari¬ 
ations  of  velocity,  temperature  and  density  for  gas 
and  particles  along  a  stagnation  streamline.  From 
Fig.  13(a),  we  can  see  that  in  the  case  of  H,  =  1.5 
and  dp  -  1  urn,  the  particle  temperature  Tp  and  gas 
temperature  T  almost  achieve  equilibrium  near  the 
stagnation  point.  However,  the  particle  and  gas 
velocities  do  not  achieve  equilibrium  near  the 
stagnation  point.  This  tendency  is  just  the 
reverse  of  the  M,  =  10  case  for  dp  =  1  urn  [see 
Fig.  “(a)).  It  is  worth  noting  that  the  variations 


ot'  the  flow  quantities  along  the  stagnation  stream- 
linos  resemble  the  known  dusty-air  solutions  behind 
a  normal  shock  wave  (Refs.  4,  b,  28,  29].  As  the 
particle  diameters  dp  increase  from  2  to  10  pm,  the 
particles  tend  to  keep  their  initial  values  of  Up 
and  Tp  as  if  the  particles  do  not  influence  the 
gas  flow. 


Shock  standoff  distance 

(a)  For  SI,  =  1.5:  the  shock  standoff  distance 
decreases  with  increasing  i  for  all  particle 
diameters  dp.  As  dp  increases,  the  effect  of 
the  freest  ream  loading  ratio  >  on  /  becomes 
weak  . 


file  gas  pressure  p  and  the  gas  and  particle 
phase  densities  .  and  . p  increase  along  the  stag¬ 
nation  streamlines  as  i  increases.  These  tenden¬ 
cies  become  weak  as  dp  increases. 

Figures  15 (a) -(d)  show  the  pressure  distribu¬ 
tions  on  the  hodv  surfaces.  As  i  increases  and 
dp  decreases,  the  pressures  on  the  body  surfaces 
become  larger.  These  effects  are  similar  to  the 
case  of  M,  =  10. 

Figures  lblal-td)  show  the  gas  and  particle 
streamlines  near  the  stagnation  regions  in  the 
shock  layers.  When  dp  -  1  .an,  the  particle 
streamlines  deflect  significantly  near  the  body 
surface  as  shown  in  Fig.  lb(n).  As  the  particle 
diameter  dp  increases  from  2  to  10  am,  the  par¬ 
ticle  deflection  decreases  and  it  proceeds  almost 
in  a  straight  line  to  collide  with  the  surface 
[Fig.  Ibid)].  This  particle  behaviour  is  similar 
to  the  case  of  M,  =  10. 


8.  CONFFUSIONS 

An  inverse  method  was  developed  to  study  dusty 
supersonic  flow  past  axisymmetric  blunt  bodies. 

The  analysis  was  based  on  the  assumptions  that 
there  are  no  collisions  between  the  individual 
particles  and  that  the  interactions  between  the 
gas  and  the  particles  are  only  through  drag  and 
heat  transfer.  The  governing  equations  for  axi¬ 
symmetric  gas-particle  flow  were  described  in  a 
cylindrical  lx ,  r) -coordinate  system.  Then  by- 
using  the  von  Mises  transformation,  the  governing 
equations  were  again  transformed  into  a  U,Y)- 
coordinate  system,  where  ,  is  the  gas  stream 
function,  and  Y  =  r.  To  facilitate  the  integration 
of  the  governing  equations  in  the  shock  layer, 
they  were  further  transformed  into  a  (-coor¬ 
dinate  system  by  using  transformation  equations 
(4.3)  and  (4.4).  The  advantages  of  the  present 
method  are:  (1)  body  shapes  are  easily  determined 
without  any  approximate  calculation:  (2)  the  com¬ 
puter  memory  required  is  not  large;  (3)  the  com¬ 
puter  calculation  time  is  very  short,  so  that  it 
is  possible  to  extend  this  method  to  a  direct 
method  by  using  an  iterative  procedure. 


(b)  For  M,  -  10:  the  shock  standoff  distance 
decreases  with  increasing  i  for  dp  -  1  urn. 
However,  when  dp  =  2,  5  and  10  i.m,  increases 
with  increasing  i,  because  of  the  decrease  of 
the  gas  density  ,  ,  in  the  shock  layer. 

Variation  of  gas  and  particle  velocities  and 
temperatures  along  a  stagnation  streamline 

(a)  For  M,  =  l.S:  when  d  =  1  nm,  the  gas  and 
particle  temperatures  T  and  Tp  almost  achieve 
equilibrium  near  the  stagnation  point  on  the 
body.  However,  the  gas  and  particle  veloci¬ 
ties  u  and  Up  do  not  achieve  equilibrium  near 
the  stagnation  point.  As  dp  increases,  the 
particles  tend  to  keep  their  initial  values 
of  lip  and  Tp. 

lb)  For  M..  =  10:  when  dp  =  1  am,  the  particle 
and  gas  velocities  u  and  Up  are  in  equili¬ 
brium  near  the  stagnation  point  on  the  body. 
However,  the  gas  and  particle  temperatures  T 
and  Tp  are  not  in  equilibrium.  This  tendency 
is  just  the  reverse  of  the  H*.  =  1.5  case. 

When  d«  =  2,  5  and  10  urn,  Up  does  not  reach 
the  value  of  u  near  the  body  surface. 


Variation  of  gas  pressure  p  and  gas  and  particle 
phase  densities  .  and  .  along  a  stagnation 
streamline  ■' 


(a)  For  H,  =  1.5:  for  all  particle  diameters  dp, 
the  gas  pressure  p  and  the  gas  and  particle 
phase  densities  t  and  v  always  increase 
along  a  stagnation  strreamline  as  <  increases. 

(b)  For  M,  =  10:  when  dp  =  1  urn,  the  gas  density 
,  decreases  slightly  along  a  stagnation 
streamline  at  first  and  then  it  increases 
near  the  body  surface.  However,  when  dp  =  2,  5 
and  10  urn,  the  density  of  the  gas  decreases 
along  a  stagnation  streamline  owing  to  the 
gas  temperature  increase  as  i  increases.  The 
density  of  the  particle  phase  ,p  and  the  gas 
pressure  p  increase  along  a  stagnation  stream- 
1  ine . 


Body-surface  pressure  p^ 


In  using  the  present  method,  the  pure-gas  flow 
fields  around  spheres  were  first  solved  numerically 
for  the  freestream  Mach  number  range  1.5  kf,.  \  10. 
These  were  found  to  be  in  very  good  agreement  with 
the  results  of  Van  Dyke  and  Gordon  [Ref.  Ib], 


For  both  kU  =  1.5  and  10,  the  body  surface 
pressure  Pb  increases  as  i  increases.  This 
tendency  is  significant  as  dp  becomes  smaller 

Behaviour  of  particle  streamlines 


Based  on  the  assumptions  that  the  Prandtl 
number  Pr  =  0.75,  the  specific  heat  of  the  gas  is 
constant,  the  specific  heat  ratio  of  the  two  phases 
is  unity,  and  the  gas  viscosity  and  drag  coefficient 
were  described  by  F.qs.  (2.20)  and  (2.15),  the  super¬ 
sonic  gas-solid  particle  flows  past  near-spherical 
bodies  were  then  solved  for  freestream  Mach  numbers 
Mo.  =  1.5  and  10,  with  freestream  loading  ratio 
a  =  0,  0.2,  0.5  and  1.0,  and  particle  diameters 
dp  =  1,  2,  5  and  10  urn,  respectively.  The  numer¬ 
ical  results  can  be  summarized  as  follows: 


For  both  M,„  =  1.5  and  10,  when  dp  =  1  nm,  the 
particle  streamlines  deflect  significantly 
near  the  body  surface.  As  the  particle 
diameter  dp  increases,  the  particle  deflection 
decreases  and  it  proceeds  almost  in  a  straight 
line  to  collide  with  the  body  surface. 

Within  the  limitations  of  the  assumptions  noted 
above,  the  present  numerical  results  provide  an 
understanding  of  dusty  supersonic  flow  past  blunt 
bodies . 
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Table  1 

Shock  Standoff  Distance  Nose  Radius  R^  and  Biuntncss  F 
of  the  Calculated  Bodies  (Spheres) 
for  Pure  Gas.  ,  *1.4 


10 

|  o 

.49 

b 

0 

.47 

4 

0. 

.38 

1 

3  1 

0 

.25 

^  ! 

-0 

.06 

1.5  i 

-0. 

,71 

Table  2 

Stagnation  Pressure  p  on  the  Blunt  Bodies 


for  Pure  Gas 


Present 

Result 

Van  Dyke 
5  Gordon 


1.01b  0.9b3  0.945  0.933  0.927 

1.007  0.957  0.940  0.928  0.923 

1.007  0.957  0.941  0.929  0.923 


Table  3 

Number  of  Solid  Particles  (Glass  Spheres) 


in  Freest  ream  Flow 


dr,  lum 


0.1 

(cm  3) 

0.2 
(cm  3) 

0.5 

(cm'3) 

1  .0 
(cm  3) 

8.98x10' 

1 .80x10® 

4.49x10® 

1 

8.98x10® 

1 . 1  2  x  1 0  7 

2.25xl07 

5.62*107 

1 .12x10® 

7.18X105 

1 .44xi06 

3.59xl06 

7 . 1  8  x 1 0® 

8.98x10** 

1  .80XJ05 

4 .49*10S 

8.98X105 

Table  4_ 

Shock  Standoff  Distance  A,  Nose  Radius  R,  and  Bluntness  B. 
of  Calculated  Bodies  for  Dusty-Cias  F 1  ows 
Mo.  =  10»  =  300  K»  =  101  *3  KPa»  =  1  cm 
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Shock  Standoff  Distance  \  Nose  Radius  R^  ami  Bluntness 
of  Calculated  Bodies  for  Dusty  lias  Mows 
M  =  l.S,  f  =  300  K,  p  =  101.3  KPa,  R  =  1  cm 
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FIG.  2  SMOCK  WAVE,  SONIC  LINE  AND  STREAMLINES  FOR  A  SPHERE 

FOR  PURE  GAS,  >  -  1.4.  »  VAN  DYKE  &  GORDON  [REF.  16]. 
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FIG.  2  -  CONTINUED 

SHOCK  WAVE,  SONIC  LINE  AND  STREAMLINES  FOR  A  SPHERE 
FOR  PURE  GAS,  y  =  1.4.  0  VAN  DYKE  5  GORDON  [REF.  16]. 
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2  -  CONTINUED 

SHOCK  WAVE,  SONIC  CINE  AND  STREAMLINES  FOR  A  SPHERE 
FOR  PURE  GAS,  >  =  1.4.  o  VAN  DYKE  S  GORDON  [REF.  16] 
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[G.  2  -  CONCLUDED 

-HOCK  WAVE,  SONIC  LINE  AND  STREAMLINES  FOR  A  SPHERE 
FOR  PURE  GAS,  Y  =  1.4.  °  VAN  DYKE  §  GORDON  [REF.  16] 


(a)  Moo  =  10,  6,  4 

FIG.  3  RELATION  BETWEEN  r2/2(x-A)  AND  (x-A)  FOR  CALCUI.ATED 
BODY  SHAPES.  o  CALCULATED  POINTS. 
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(b)  M®  =3,2,  1.5 


FIG.  3  -  CONCLUDED 

RELATION  BETWEEN  r  /(2(x-A)  AND  (x-A)  FOR  CALCULATED 
BODY  SHAPES.  °  CALCULATED  POINTS. 


CONTINUED 

PRESSURE  DISTRIBUTION  ON  THE  BODY  (SPHERE)  SURFACE 
FOR  PURE  GAS,  y  =  1.4.  0  VAN  DYKE  §  GORDON  [REF.  16] 
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FI Ci.  4  -  CONTINUED 

PRESSURE  DISTRIBUTION  IN  THE  BODY  (SPHERE)  SURFACE 
FOR  PURE  GAS,  >  =  1.4.  o  VAN  DYKE  5  GORDON  (REF.  ID], 


FIG.  5  ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES  FOR 
M  =  10,  T  =  300  K,  p  =  101.3  KPa  AND  R=  =  1  CM 
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5  -  CONTINUED 

ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES  FOR 
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FIG. 
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5  -  CONTINUED 

ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES  FOR 
=  10»  =  300  K,  =  101.3  KPa  AND  Rs  =  1  CM. 
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(d)  dp  =  10 fim 


-  CONCLUDED 

ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES  FOR 
Ho  =  ln<  =  -™0  K,  Pm  =  101.3  KPa  AND  Rs  =  1  CM 
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FIG.  ?  -  CONCLUDED  7 

RELATIONS  BETWEEN  r“/2(x-A)  AND  (x-A)  FOR  CALCULATED 
BODIES  FOR  Mx,  =  10,  =  300  K,  =  101.3  KPa  AND 


U,  Up  ,T ,  Tp 
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(a)  dp=  l/*m 


FIG.  7  VARIATION  OF  GAS  AND  PARTICLE  VELOCITY  AND  TEMPERATURE 
WITH  LOADING  RATIO  a  ALONG  STAGNATION  STREAMLINES  FOR 
H„  =  10,  TttJ  =  300  K,  P<„  =  101.3  KPa  AND  Rs  =  1  CM. 
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(c)  dp  =  5/i.m 


FIG.  7  -  CONTINUED 

VARIATION  OF  GAS  AND  PARTICLE  VELOCITY  AND  TEMPERATU,,.. 
WITH  LOADING  RATIO  <x  ALONG  STAGNATION  STREAMLINES  FOR 
H,,  =  10,  T„,  =  300  K,  p[X>  =  1 0 1  . 3  KPa  AND  Rs  =  1  CM. 
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(c)  dp  =  5^m 


FIH.  8  -  CONTINUED 

VARIATION  OF  (IAS  AND  PARTICI.F.  DF.NSITlIiS,  AND  GAS  PRKSSURF 
ALONG  STAGNATION  STRLAMLINF.S  FOR  Mu  =  10,  f  =  300  K, 
p  =  101.3  KPa  AND  Rs  =  1  CM. 
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8  -  CONCLUDED 

VARIATION  OF  GAS  AND  PARTICLE  DENSITIES,  AND  GAS  PRESSURE 
ALONG  STAGNATION  STREAMLINES  FOR  =  10,  Tro  =  300  K, 
p  =  101.3  KPa  AND  Rs  =  1  CM. 
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ON  BODY  SURFACE  FOR  M  =10 
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CONTINUED 

PRESSURE  DISTRIBUTIONS  ON  BODY  SURFACE  FOR  M 
p„,  =  101.3  KPa  AND  Rs  =  1  CM. 
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FIG.  9  -  CONTINUED 

PRESSURE  DISTRIBUTIONS  ON  BODY  SURFACE  FOR  MOT  =  10,  T 
p_  =  101.3  KPa  AND  Rs  =  1  CM. 


=  300  K, 
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(a)  dp  =  1/im,  a=  0.2 


FIG.  10  GAS  AND  PARTICLE  STREAMLINES  AROUND  BLUNT  BODY  FOR 
Mco  =  10*  T<„  =  300  Pa,  =  101.3  KPa,  Rs  =  1  CM. 

-  GAS  STREAMLINE,  -  PARTICLE  STREAMLINE, 

IMPACT  ANGLE  6. 
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(b)  dp  =  2/xm,  a  =  0.5 


FIG.  10  -  CONTINUED 

GAS  AND  PARTICLE  STREAMLINES  AROUND  BLUNT  BODY  FOR 
M  =  10,  f  =  300  K,  p  =  101.3  KPa,  =  1  CM. 

-  GAS  STREAMLINE,  -  PARTICLE  STREAMLINE, 

IMPACT  ANGLE  A. 
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(c)  dp  =  5  /irr\ 


1  -  CONTINUED 

ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES 
M~  =  1’5*  Too  =  300  K>  P„,  =  101.3  KPa ,  Rs  =  1  CM. 


/ 


_  X 
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FIG.  11  -  CONCLUDED 

ASSUMED  SHOCK  WAVE  AND  CALCULATED  BODY  SHAPES. 
Moo  =  1-5,  =  300  K,  =  101.3  KPa ,  Rs  =  1  O 
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(b)  dp  =  2/im 


FIG.  12  -  CONTINUED 

RELATION  BETWEEN  r  /2(x-A)  AND  (x-A)  FOR  CALCULATED  BODIES. 
Mcc  =  1-5*  -  3°°  =  101.3  KPa,  Rs  =  1  CM. 
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FIG.  13  VARIATION  OF  GAS  AND  PARTICLE  VELOCITIES  AND  TEMPERATURE 
ALONG  STAGNATION  STREAMLINES.  M  =  1.5,  f  =  300  K, 
pOT  =  101.3  KPa,  Rs  =  1  CM. 
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FIG.  13  -  CONTINUED 

VARIATION  OF  GAS  AND  PARTICLE  VELOCITIES  AND  TEMPERATURE 
ALONG  STAGNATION_STREAMLINES .  M^  =  1.5,  T^  =  300  K, 

=  101.3  KPa,  Rs  =  1  CM. 


FIG.  13  -  CONCLUDED 

VARIATION  OF  GAS  AND  PARTICLE  VELOCITIES  AND  TEMPERATURE 
ALONG  STAGNATION  STREAMLINES.  M  =15  f  =  300  K 
=  101.3  KPa,  Rs  =  1  CM. 


(a)  dp  =  j  /xm 

FIG.  14  VARIATION  OF  GAS  AND  PARTICLE  PHASE  DENSITIES  AND  GAS  PRESSURE 
ALONG  STAGNATION  STREAMLINES.  M  =  1.5,  T  =  300  K,  p  =  101.3 
KPa,  Rs  =  1  CM. 
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FIG.  14  -  CONTINUED 

VARIATION  OF  GAS  AND  PARTICLE  PHASE  DENSITIES  AND  GAS_PRESSURE 
ALONG  STAGNATION  STREAMLINES.  M  =  1.5,  T  =  300  K,  p  =  101.3 
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FIG.  14  -  CONCLUDED 

VARIATION  OF  GAS  AND  PARTICLE  PHASE 
ALONG  STAGNATION  STREAMLINES.  M  = 
KPa,  Rs  =  1  CM. 


DENSITIES  AND  GASPRESSURE 
1.5,  T^  =  300  K,  =  101.3 
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(a)  dp  =  l/tm,  a  =  0.2 

=IG.  16  GAS  AND  PARTICLE  STREAMLINES  AROUND  BLUNT  BODY. 

M  =  1.5,  f  =  300  K,  p  =  101.3  KPa,  Rs  =  1  CM. 
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-  GAS  STREAMLINE,  -  PARTICLE  STREAMLINE. 
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FIG.  16  -  CONTINUED 

GAS  AND  PARTICLE  STREAMLINES  AROUND  BLUNT  BODY. 

M  =  1.5,  T  =  300  K,  p  =  101.3  KPa,  Rs  =  1  CM. 
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FIG.  16  -  CONCLUDED 

GAS  AND  PARTICLE  STREAMLINES  AROUND  BLUNT  BODY 
M=o  =  15*  foo  =  300  K,  P„  =  101.3  KPa,  Rs  =  1  CM. 
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